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1  INTRODUCTION 


-  -  The  aim  of  one  proposed  investigation  was  to  generate  a  set  of  theoretical 
and  experimental  unsteady  aerodynamic  data  for  an  advanced  type  -airfoil.  I  his  set 
of  data  should  serve  as  a  standard  for  the  comparison  and  the  evaluation  :f  com¬ 
putational  methods  for  two-dimensional  unsteady  transonic  flow.  For  this  turpose 
it  was  intended  to  complement  the  results  of  the  wind  tunnel  tests  conducted  at 
NLR  on  the  supercritical  NLR  ”301  airfoil  with  theoretical  results,  to  he  computed 
with  the  method  of  Dr.  R.J.  Magnus  and  Dr.  H.  Yoshihara  *  ot'"j5ners.l  Dynamics, 

San  Diego.  The  calculations,  to  be  performed  'ey  General  Dynamics,  v  ere  supposed 
to  be  sponsored  'by  the  office  of  Naval  Research 
*  The  statement  of  work  of  the  intended  investigation  was: 
mu  NLR  provides  the  contour  da' a  of  the  "shock-free"  NLR  7301  airfoil  sect:  on,  to¬ 
gether  with  the  hodo graph  solution  for  rhe  design  condition.  • 

-  NLR  provides  the  measured  steady  and  unsteady  pressure  distributions  ir.  three 
characi  eristic  flow  concisions  :  subsonic  flow,  transonic  flow  w:.tn  a  well  de¬ 
veloped  supercritical  region  terminated  by  a  shock  wave  and  the  "shock- 're e" 
design  condition. 

-  General  Dynamics  performs  computations  with  their  "exact"  method  to  obtain  the¬ 
oretical  results  for  the  three  mentioned  flow  conditions. 

-  General  Dynamics  and  NLR  will  put  together  the  corresponding  results  in  a  com¬ 
mon  report,  including  a  critical  evaluation. 

Unfortunately,  due  to  lack  of  sufficient  funds  to  perform  the  final  sot  of 
computations  the  proposed  program  could  not  be  completed 

This  report  gives  a  summary  of  the  work  performed  so  far. 


2  EXPERIMENTS 


At  NLR  an  extensive  wind  tunnel  investigation  has  been  conducted  to  explore 
the  unsteady  aerodynamic  characteristic  of  sc— called  supercritical  airfoils.  In 

this  framework  detailed  steady  an  unsteady  pressure  distributions  have  been  _ 

measured  on  the  oscillating  16.5  percent  thick  NLR  7301  airfoil  (  for  contour 
data  see  table  l),  which  was  designed  with  the  hoaograph  method  of  Boerstoel 


(  Ref 


i  \ 


Codes 

l/or 


Jlsi 


bpooiai 


A 


*  J 


(*)  At  present  Dr.  Yosnihara  is  employed  by  the  Boeing  Co 


□  □ 
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Emphasis  was  put  on  three  mean  flow  conditions,  which  can  be  characterized  as 
follows  : 

A  :  subsonic  flow 

B  :  transonic  flow  with  a  well  developed  supersonic  region  on  the  upper  surface 
terminated  by  a  relatively  strong  shock  wave 
C  :  the  "shock-free"  design  condition. 

According  to  the  statement  of  -work,  the  contour  data  and  the  test  results 
were  made  available  to  General  Dynamics.  Further,  independently  of  the  present 
investigation,  the  results  of  the  experiments,  including  a  thorough  analysis, 
have  been  presented  during1  the  AGARD  conference  on  ri  narto  -in  spiwm- 

ted  and  transonic  flow" ,  Lisbon,  April  1977-  For  convenience  of  the  reader,  a  re¬ 
print  of  the  AGARD  paper  is  added  to  this  report  as  Appendix  A. 

Meanwhile,  the  NLR  7301  airfoil  has  been  selected  one  of  the  test  cases 
for  the  AGARD  activity  "Standard  Aeroelastic  Configurations".  Fcr  that  purpose 
all  relevant  data  (geometry  description,  aerodynamic  conditions  and  tabulated  ex- 
perimenlal  values)  have  been  gathered  in  tabular  form  in  reference  2. 


3  COMPUTATIONS 

As  stated  earlier,  it  was  intended  to  complement  the  results  of  the  experi¬ 
ments  with  results  obtained  with  the  computational  method  developed  by  Magnus  and 
Yoshihara  (Refs.  3.  4).  Their  method,  which  can  be  considered  as  one  of  the  most 
advanced  methods  at  prese  vailable  for  unsteady  flow  computations,  solves  the 
Euler  equations  without  fv.  1  assumptions  concerning  frequency  or  amplitude?  of 
oscillation.  The  boundary  conditions  are  imposed  along  a  contour  coincident  with 
the  mean  position  of  the  airfoil.  The  computations  start  with  the  steady  flow 
solution  and  are  continued  until  a  complete  cycle  of  the  periodic  flow  is  obtai¬ 
ned.  Til  is  is  a  time-consuming  process  and  therefore  the  method  is  not  suited  for 
routine  use,  but  rather  to  generate  some  solutions  that  might  reveal  the  nature 
of  the  flow  or  serve  as  test  cases  for  more  approximative  but  faster  methods 

Dr.  Magnus  performed  a  set  of  computations  for  the  NLR  7301  airfoil  and  ga¬ 
thered  the  results  in  report  CASD/LVP  78-013  entitled  "  Some  numerical  solutions 
of  inviscid,  unsteady  transonic  flow  over  the  NLR  7301  airfoil"  (Ref.  5)*  His 
results  showed  many  of  the  features  of  the  experiments,  but  that  there  were 
differences  due  to  the  effect  of  the  tunnelwalls,  the  presence  of  the  boudary 
layer  and  the  fact  that  in  the  calculated  results  insufficient  expansion  of  the 
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flew  along  the  upper  part  of  the  nos°  could  be  obtained  (Fig.  l)  .  Further  it  was 

experienced  that  the  final  result  was  influenced  considerably  when  m  the  compu¬ 
tations  the  trailing  edge  was  chopped. 

In  view  of  these  results  it  was  decided  that  the  present  study  could  only  be 
finished  in  a  fruitful  way  if  the  following  additional  calculations  were  perfor¬ 
med  : 


condition 

M 

geometric 

incidence 

ao 

incidence 
corrected 
for  wall 
interference 

ac 

reduced  frequency 

wc 

*  ”  2U 

A1 

0.5 

(0.35°) 

0.40° 

0.263 

(3.00°) 

o 

4ith  tail 

B1 

0.7 

2.00 

0.192 

Cl 

0.721 

- 

-0.19° 

0.1S1 

condition 

M 

a 

a 

k 

0 

c 

B 

0.7 

(3.00°) 

2.00° 

0.192 

tail  chopped 

c 

0.721 

-  0.19° 

0.181 

at  x  =  1.00 

The  amplitude  of  oscillation  around  the  pitch  axis  at  40  percent  chord  amounts 


Essential  is  that  the  computations  were  to  be  performed  for  mean  incidences 
ac  in  which  the  effect  of  the  tunnel  walls  is  taken  into  account.  In  this  way  it 
is  assured  that  in  the  calculations osci llatory  perturbations  around  the  same 
mean  steady  flow  field  are  considered  as  occurred  in  the  experiments. 

As  far  as  the  "shock— free"  flow  condition  is  concerned,  preference  is  given  to 
the  theoretical  design  condition,  which  is  found  at  M  =  0.721  and  a  =  -  .19C 
Unfortunately,  the  computations  for  the  cases  summari zed  m  the  foregoing 
table  cculd  not  be  performed,  since  no  funds  could  be  made  available  to  Dr. 
Magnus  to  make  the  additional  computer  runs. 
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4  CONCLUDING  REMARKS 

Due  to  the  fact  that  the  final  sex  of  computations  cannot  be  performed,  al¬ 
so  the  remaining  tasks  for  NLR,  namely  a  comparison  between  theory  and  experiment 
and  a  critical  evaluation  of  tne  results  in  collaboration  with  General  Dynamics 
cannot  be  completed.  The  APOSR  funds  which  were  made  available  for  this  part  of 
the  NL\  work  will  be  returned  to  APOSR. 
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TABLE  1 

Contour  data  of  airfoil  NLR  7301 


UPPER  PART 


lower  par  t 


X/C  Z/C  x/c  z/c 


x/c  z/c  x/c  z/c 


0.300001? 
O.O0O289S 
0. 0008873 
0.1)008*61 
0.0011758 
ft./>1l4602 
0.0017440 
3.  on?3  1  7 s 
0.602298“ 
ft. 1028701 
ft . Oft?*  7 1  * 
0.ftft4??04 
ft. 3086121 

3.3  349?  39 
3 . 304 l 44? 
ft. ft39* 70S 
0. ft  10 78. 7ft 

ft. 19447 

3.01311** 
3.3)4>9)P 
ft.ft 1830*9 

ft .31443*4 

3.3178313 
3.31 *874  ft 
0.3194ft*? 
ft.ft?ft>?44 
3.3?l 1 ?  7  3 
ft.n??ft*ft? 
ft.  3^40880 
3.3?88**7 
ft .  0?482 7 4 
ft  .  3?7490« 
ft.4?7 4023 
0.3274*3“ 
0.3279647 
1 .  " 28 1 ‘ n  • 
ft. 0208*19 
0.3791 330 
0.0294602 
3.0700488 
0.0304510 
0.311049? 
3.0314*79 
0.3331274 
0.0378378 
0.04284*3 
0.0472590 
0.0519708 
3  •  0  566*?  4 
■V*  !  10*4 
1.3441071 
ft  .  *  *  rt  A  j  O  a 
1.-->74400-> 

3.3*4**?* 

'  .  .ft  *  *  ■  0  9  7 
ft.H?4S8* 
3  ft  3*.  7  1-  7 
ft.  1  0?*7<J  1 
3.  1 

ft.  117-44) 
,.  \  334Q11 
ft. i  79*  3  7 « 
ft.  1  '’l’^ 

0  •  i  4  4  (  4  |  w 

3,  ) 8 1 1 4  9  J 

3.14><)?  .7 
3. 1 7391 0ft 
ft. 1 78741 s 
ft . 1 *131 34 
0.  J  90  ?  7  5  7 
ft,  1  9*4  3  7<« 
0.?Q39991 
1,  ?1 4 1A/H 

0.711 7??4 

3  ,  > \  7497  ft 
0,  >73471  4 
3. *791459 
o .  7  352?0S 

3.2417117 
0.24*?030 
). 2844944 

’  .  7411*5* 
3.  ?4*165l 
0. *78844 J 
1, ?*27?34 
3,?*990?9 
rt.  P97T7SJ 
n.  1086474 
0.11 18?3ft 
3.  1213928 

3.1303*47 

3.11*7013 
0.1473884 
0.  1840399 
0.1881414 
0.  IF*?  7.14 
0.ft*34081 
0. J925371 
0, *91 *448 
0.411 1518 
0.410489? 


-.0.304  14? 
*.0fi5?l91 
•.007476* 
• .0098948 
•.01 122) 7 
..0178337 

•  .0 1 369  ft* 

•  .31*04  1  1 

♦.C16087A 
*.riTori4 
..Ul 79*7“ 
2?05? * 
.  7  ? 

».<■?*  1  7ft  > 

.. ’10831? 
..0 1266  7? 
•.0346801 
..C16  1896 

• .037*449 
•.018 | 4 7s 
..0404096 
.. 041403s 

•  . 0*2681* 
..**.161*7 
.  .  .  *  1**  ? 

'•.4  7Q1  > 
..  .673ft.ft 

»  .  0*  7ss-*  1 

».  34**03^ 

«  .  0.01749 

•  .  0.0904* 

•t0. 00651 
.  .09AS-. 

• . 050 1 55  ft 
*.3432416 

•  .  ;;  4  0 . 7  J  6 
..rs37.»\ 

•  .051005? 
*.051191“ 
*.051 3840 

• . 0e 1 6494 
•.081970? 

•  .  -ft  4  7  6  1 4  1 
..08*6577 

•  .  J  6  6  4  7  ?  3 

•  .04*106* 
.  .  0«94*2< 

•  .0609237 
* .  '6  2;  5  74 

•  .06328' * 
..064  13  ftft 
-.064.949 
..0*6  17  0s 
• .067135* 

.  . 347*537 
.  .06*5241 

•  .  0  6  8  7  1  S  9 

•  .0  702183 

•  ..071  14*4 

-  .0723984 

.  .  30979 

•  .0  7.0  ’  3* 
..178. .Ji* 

•  .  0  7*959 7 

..ft»6*^7t 

..«»?*??? 

•  .  07*  7?7Ci 

..07919] * 

-  .  079*469 
* .  0*381 3 3 

-  .  ''■3945* 

•  .0*1 *343 
*.r *1 *?4T 
..0*2231* 
..0*76671 
..0*13779 

.  . c*  1*449 

•  .0*38* l “ 

•  .0**217  7 

..0*461 35 

•  .3*4  963* 
»  .0*8  3068 
.  .  0*44577 
»  .0*49889 
t/PAJQja 
.  .0*68787 
.  .  0*4167? 

•  .0*71 229 
..0*715*4 
•.0*75690 

•  .  0*77771 
•.0079*41 
•.0**0910 
*.O*0?O7i 
..0**2986 
..00*358* 
..00*3067 
..00*3033 

•  .0**3*75 
..00027*9 

•  •  3** 1 778 


ft.  *29746  7 

..l>“ft*  1* 

.  /  0  J  V  1  2 

-  .  .  ft./  -  1  6  2 

9 . -s--s  74 

-  0728)26 

3.4  1*43*4 

• . '*7***7 

-  .0)02  J ‘ift 

-  .  UCS62  J  ’ 

9 . -008 1 8  7 

-  Q717100 

0.4474*29 

. . 0“7s**^ 

-  .  J  ft  >0  2  t 

.  .  -  7-6000 

-  07H/08 

3.4447*11 

. . '“7*41* 

- . ft  1C  1 ib* 

- . 00-0- 3- 

9  .  -844 ) )? 

-  0700207 

0 . *48  »  1 9  > 

- .  47?)7* 

.  S-'O  —  j , 

-.9  1  J2297' 

9.—*) J97| 

-  90*1612 

ft . »)/». 1 1 

. . - 9*99 

>  .3.1  1  7  *-> 

:  1 121-* 

..8001  7  0S 

-  0603)70 

0.*“04ft3ft 

• . '•‘J?)- 

.  12»9 > 

9.8007021 

-00781 38 

.  .  -  OS.  17  » 

I.V9 iftftl* 

-.91  2*7  •»  J 

1. 5l))-*7 

-  .  06O636  7 

3. ,943670 

..'“•11?: 

-.0.1 77*- 

- . . 1  3)2-9 

1  • 82022-h 

- . 0046*3* 

ft  ,4ft7s1 1  > 

..••“4*1*7 

•  •  J  9  1  l*  *  f  - 

■.Jli'HJ 

J  .  8?  7-883 

-  .  06-6603 

ft  ,  8  0  O  1  ft  4  ft 

.  .  w46  >s  . 

..ft-Vlftift 

- . J 1 --6. fl 

-  .4  )8 1272 

-.06  )8-  l  - 

1  •  8  l  8  7  1  9  1 

:s 

-  .  JV- 

.  .  9  1  s  1  '9  * 

9.5-317  — 

ft02J) *1 

.1.472761? 

. .  c*  AS  )  1 

. .  ft  ft  J  1  >*> 

- . J i s*  .  - 

J.  44  l  44  7  9 

- .0009*  '  3 

0 .8  11797  4 

. .  “** ' 

-  -  0  3  -  1  ) 

-  .  1’  J  -  «  1  s  - 

9.58-4-** 

0&U423- 

0.51*  111  1 

„  .  ft  )-1'  1  3  > 

-.01 *76  )- 

.8838107 

-.94*0  )  7  j 

'  . 5*  14*“  » 

. .  0  '■'ye  747 

- .Or  3*7 *J 

9.8'?5«3- 

-. ftS  780  )7 

ft. 8575 '7 > 

.  .  0  ft  ft  **  i  i  J 

-  .  9  7  1  *  ?  2  * 

. 482 1 8sb 

)884ftl- 

ft  ,  5*. }  I?-; 

. .  a?  I  .  t  . 

- .  ft^esr-b 

•  .  .  t  i  0  -  5  » 

.8*I62*C 

- . >4— 1 -83 

C .84*447  » 

. .  ’  0  **  »  > 

-  .0  j»  «**? 

-  .  9  '-9S-- 

.59 | 1 U82 

-.  )42  3*  3- 

1.4749*9  . 

•  .  •  1  >  7  :  • 

.  .  i-  0 

- .  -2s :  i  i' 

.6104822 

- . 14O4802 

ft ,4J  7?>s  « 

.  .  ■**  >>-  ’ 

.ft  >1  .  - 

-  .  -  74  . ss 

.61*72-  ) 

-  .  Uo'no 

*  .  5*ss4w  . 

. .  «-»;'! 

.  .  j  ft  -  •'*,  1 

-  .  9  2  ’ ;  i  -  ft 

. 8208682 

- .  )-0  *60 1 

'  .  -.J)«4k  » 

•  .  70  7  1  >  . 

.  .  C  1  ft  "  2  i  ' 

-  .  .  2  <*  v  2  *  ' 

.8176224 

-  .  1-8I 6*  — 

1.5 >974*  ; 

..  ‘  T  O  *  7  9  7 

-.011  '7*0 

-  .  9  c  -  *  - “ 

.8-6  1 706 

)-  >  3  )66 

'  .  4  ft  .  '  •.  J  - 

.  .  r  “  1  '  '  A 

.  .0  1  2  7*»4J 

-  •  <•*-*-- 

.  8sS  J8  )4 

-  ■  1-  1  -88* 

•' ■  .  4  ft  *  7  ,  -  | 

•  • 

.  l  J  * . :  * 

-  .  J  ,  -  s  '  - 

.05372«3 

- .  J*6 1 -  * 
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A/3532 


Date 
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Dear  Capt.  Wilkins, 


Thank  you  for  your  letter  of  June  19th,  dealing  with  the  AROSE 
"Unsteady  airloads  or.  a  sinusoidally  oscillating  supercritical 
foil",  which  expired  cn  Dec.  31st.,  1978. 


Grant 

air- 


As  you  know,  she  aim  cf  the  proposed  investigation  was  to  generate  a 
set  of  unsteady  aerodynamic  data  for  an  advanced  type  of  airfoil, 
noth  theoretical  and  experimental,  which  may  serve  as  a  standard  for 
comparison  and  evalvrm ion  of  calculation  methods  for  two-dimensional 
unsteady  transonic  flow.  For  this  purpose,  it  was  intended  to  comple¬ 
ment  the  results  of  the  wind  tunnel  tests  conducted  at  the  NLR  cn  the 
lo.p  percent  thick  NLR  7301  airfoil  with  calculated  results  obtained 
with  the  method  of  Dr.  Richard  J.  Magnus  and  Dr.  H.  Joshihara  of 
Her. oral  dynamic.-.  Gan  Diego.  In  collaboration  with  General  Dynamics 
a  comparison  and  critical  evaluation  should  be  made  of  both  theory 
and  experiment . 

The  calculations  to  be  performed  by  General  Dynamics  were  supposed 
to  be  sponsored  by  the  Office  of  Naval  Research. 

MLR  provided  the  contour  data  of  the  airfoil  and  the  relevant  test 
data  to  General  Dynamics.  Further,  apart  from  the  present  investiga¬ 
tion,  an  analysis  of  the  test  results  was  presented  during  the  AC-ARD 
Conference  on  "Unsteady  airloads  in  separated  and  transonic  flow", 
Lisbon,  April  1r'77  (see  Agora  CP  226).  Moreover,  details  of  the  ex¬ 
periments  and  its  results  have  been  published  in  NLR  report  TR  77090  U 
"Investigations  of  the  transonic  flow  past  oscillating  airfoils". 

Meanwhile,  Dr.  Magnus  performed  a  first  set  of  calculations  on  the 
NLR  7301  airfoil  and  sent  us  some  results  with  his  letter  of  12  Jan. 

./.  1978  (see  enclosure  A).  This  set  of  results  has  been  published  in 
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unsteady,  transonic  flews  the  NLR  73  '1  airfoil"  'tv  R  .  Magnus 

(1078).  The  Tester  of  P r .  Magnus  (enclosure  Aj  gives  a  recur r.e  of  the 
unforeseen  problems  ho  encountered  during  the  computations  and  his 
financial  problems  in  continuing  the  work . 

In  our  letter  of  April  11  to  (enclosure  B),  v»  rave  some  comment s 
cn  the  results  ccmpu*  cu  so  far  and  ir.  Heated  that  :n  order  to  continue 
our  part  of  trie  :  over’. igaticns  in  a  fruitful  vav ,  some  additional  ral  - 
culatior.s  verc*  required. 


These  calculations  van---  discussed  in  detail  between  Pr.  Magnus ,  Pr . 
Tijdexan,  and  Pr.  Par.  for  i  lav  is  of  IIASA  during  th.e  workshop  on  transonic 
aerodynamics  for  aercelastin  applications,  ( Columbus ,  Ohio,  Ncv .  1978). 
Agreement  vas  r^ach-ei  on  tire  specific  calculations  still  to  re  per¬ 


formed.  Ar.  inter-: 


brought  up  by  Pr.  Magnus,  van  the  sensi¬ 


tivity  of  the  results  to  the  cut-off  of  the  trailing  -  i re  and  it  vas 
decided  also  to  add  seme  calculations  to  reveal  this  effect. 


_  ro  gram  vas  e->nf  irm.-'d  in  cur  le*  ter  ~f  11  Pan.  If:  (see  en¬ 
closure  C)  and,  in  cri--r  to  avoid  any  misunderstanding,  a  set  of  the 
airfoil  contour  data  vas  enclosed  again. 


r rem  the  letter  ve  received  from  Pr.  Davis  (see  enclosure  Pi  ve  learned 
that  no  funds  could  be  provided  to  perform  the  final  set  of  calculations 
by  Pr .  Magnus . 


If,  vi ev  cf  this  devils rment  it  is  not  possible  for  us  to  finish  the 
work  as  originally  planned.  For  this  reason  ve  ask  you  formal  permission 
tc  stop  the  work  under  the  Grant.  If  you  agree  ve  will  write  a  brief 
scientific  re per?  and  send  it  to  major  low ell  at  HOARD  together  vi*h 
a  final  invent. ion  report  and  a  final  fiscal  report  as  noon  as  rossible . 

A  check  for  the  uns;  or.t  grant  funds  will  be  sent  directly  to  your  office. 


In  our  opinion  it  is  .a  pity  that  the  work 
cussed  this  vith  Dr.  Panford  Davis  (KASA) 
and  they  feel  the  same.  During  a  visit  of 
discussed  the  following  asT.,e-_ls: 


cannot  be  completed.  Me  dis¬ 
and  Mr.  Charles  F.  Coe  (NASA) 
Mr.  Coe  at  NLR  last  June  we 


-  Dr.  Chyu  of  NASA  Amos  is  able  to  perform  similar  computations  for 
the  NLR  7301  airfoil  as  originally  planned  by  Dr.  Magnus. 

The  computer  program  of  Dr.  Chyu ,  however,  is  much  faster. 

-  NASA.  Ames  has  performed  an  experimental  program  on  the  oscillating 
NLR  7301  airfoil  in  their  11  ft.  transonic  winitunnel  putting 
emphasis  on  the  effect  of  Reynolds  number,  amplitude  of  oscilla¬ 
ting  and  mode  of  vibration. 
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i:-: 

1.  INTRODUCTION 

Nowadays  there  is  a  considerably  interest.  m  methods  to  pr>M:ei.  the  unsteady  air.  .ads  air¬ 
foils  and  wings  oscillating  in  transonic  flow,  '..specially  in  connection  with  the  current  interest  in 
the  so-called  supercritical  wir.g  concept.  H.wwer,  in  contrast  with  trie  stea  y  flow  cane,  experimental 
data,  that  are  sufficiently  detamed  to  verify  fundamental  theoretical  assumptions  or  to  confirm  tne 
validity  of  calculated  results  are  very  scarce  and  thus  a  definite  need  exists. 

For  this  reason  recently  at  NLR  an  exploratory  wind  tunnel  investigation  nas  seer,  performed  or. 
a -model  of  an  oscillating  supercritical  airfoil,  of  which  the  geometry'  iias  Leer  generated  with  trie 
hodograph  method  of  Boerstoel  (Refs.  1,  2).  While  the  airfoil  was  oscillating  in  pitch  about  an  axis 
at  liO  per  cent  of  the  chord  detailed  pressure  distributions  were  determined.  Ii  addition  time  his¬ 
tories  of  shock  wave  motions  were  recorded. 

The  aim  of  the  present  paper  is  to  illustrate  some  typical  high  sui  sonic  ar.u  transonic  effect:  as 
observed  in  the  experiments.  After  a  brief  description  cl  the  test  set  up,  an  analysis  is  given  of  the 
pressure  distributions  and  the  resulting  unsteady  airloads  as  measured  for  some  characteristic  flow 
conditions.  Further  attention  is  paid  to  the  periodical  motions  of  the  shock  wave  and  finally  it  is 
tried  to  assess  wnat  can  be  expected  from  the  new  generation  of  calculation  methods  for  unsteady 
transonic  flow  iFor  details  about  the  various  theoretical  methods  reference  is  made  to  the  otr.er 
papers  presented  during  this  meeting). 

2.  MODEL  AND  TEST  SET  UP 

2.1  Model  and  excitation  system 

The  airfoil  under  consideration,  the  NLR  7201,  was  designed  for  "shock-free"  flow  under  pres¬ 
cribed  conditions  (Fig.  1)  and  was  tested  extensively  -r.  steady  flow  by  Hohr.e  and  Zwaaneveid  (Refs,  i, 
1).  For  the  purpose  of  the  present  unsteady  experiments  a  new  model  has  beer,  built,  which  could  per¬ 
form  pitching  oscillations  about  ar.  axis  at  -0  p°r  cent  ;f  the  chord.  This  model,  maae  of  Dural,  has 
a  chord  length  of  15  cm  and  spans  heritentaiiy  the  test  section  of  the  NLR  Pilot  tunnel.  The  pitching 
motion  is  generated  by  means  of  a  hydraulic  actuator  (For  a  detailed  description  of  the  hydraulic 
system  and  the  model  suspension  reference  is  made  to  Poestkcke  (Ref.  5)),  To  keep  the  suspension  as 
simple  as  possible  the  model  is  excited  at  one  side,  whil  =  the  opposite  side  is  supported  by  a  bearing 
just  outside  the  tunnel  wall  (Fig.  2).  To  avoid  acomp.icated  sealing  between  model  and  window,  the 
window  closest  to  the  actuator  is  attached  to  the  model  and  fellows  its  notion.  In  addition  it  results 
in  a  clear  view  or.  the  model  surface  for  the  optical  flow  studies. 

Both  the  upper  and  lower  surface  of  the  model  are  provided  with  20  pressure  orifices  (rig.  31, 
connected  with  two  scanning  valves  outside  the  wind  tunnel  via  pressure  tubes.  Ir.  addition  11  minia¬ 
ture  Kulite  transducers  are  built  ir..  This  number,  which  is  larger  than  necessary  for  the  dynamic 
calibration  of  the  pressure  tubes,  was  chosen  to  create  the  possibility  to  arrest  the  actual  time 
histories  (including  the  higher  harmonics)  of  the  chord-wise  pressure  distribution  along  the  upper 
surface. 

To  determine  the  motion  of  the  model  use  is  made  of  6  accelerometers,  located  ir.  three  spanwise 
sections.  The  mean  incidence  is  controlled  by  the  hydraulic  system. 

2.2  Optical  flow  studies 

The  periodical  shock  wave  motions  or.  the  oscillating  model  were  determined  from  a  series  of  sub¬ 
sequent  shadowgraph  pictures.  These  pictures  were  tak"r.  using  a  stroboscopic  light  source,  trigger  d 
by  an  electrical  signal  from  a  displacement  rick  up.  By  means  of  ar.  adjustable  phase  shift  in  the 
electric  circuit  between  the  accelerometer  and  the  ligr.t  source  the  oscillating  model  with  its 
instantaneous  shock  pattern  cou-d  be  photograpned  in  every  position  desired. 

2.3  Wind  tunnel 

The  experiments  were  performed  ir.  the  Pilot  tunnel  of  i.'LR,  which  is  an  atmospheric  closed  circui 
tunnel  for  Mach  numbers  up  to  1.  Upper  and  lower  surface  of  the  test  section  'height:  55  cm;  width: 

-2  cm/  are  fitted  witr.  longitudinal  slotted  walls.  The  oper.  area  ratio  of  the  walls  is  D.l  and  the 
plenum  chambers  of  fleer  and  bettor,  are  not  connected.  Further  details  of  the  Pilot  tunnel  can  be 
found  ir.  reference  r . 


3.  ANALYSIS  IF  RESULTS 
3.1  Introductory  remarks 
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ponding  quasi-steady  airload 
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ami:  characteristics  of  the  NLP  7301  airfcil  will  be  discussed  using 
fferer.t  flew  conditions,  which  can  be  characterised  as  follows  (see 


shock  wave 


the  dynamic  effects  or.  the  unsteady  airloads, 
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where  Aa  is  the  amplitude  of  the  pttcr.ir.g  csci llaticr.  ir.  radians,  p  the  treasure  variation  and  p  the 
dynamic  pressure.  The  correct  end  mg  :caoi-cteady  coefficient  can  le  derived  from  steady  tests  at  two 
differ®.-.*  ;r.c;  ier.ee 
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•ingle  of  rit  t  aex  *i:.u  :  quas  i-.-t  *"i-!v  result;;  :°r  iv*vi  from  f  at  an?  rnovr.  in  f  i  gu  re  r* .  In  order  t  ? 

facilitate  tho  comparison  between  upper  ana  i.*v^r  rurfaee  *:.«  quas i-steady  pressures  at  the  urp^r 
surface  are  plotduu  with  a  reversed  sign  v  F  i  g .  '  .II),  I:.e  agreement  bet  veer,  tne  measured  quasi-steauy 
pressure?  and  th°  predict:  or.  of  rbir.  airfoil  theory  is  reasonable.  The  largest  deviations  show  up 
over  the  rear  part  of  *  he  airfoil,  where  the  measured  data  are  below  tne  calculated  curve. 

Near  the  leading  edge  the  measured  pressures  on  tne  urper  surface  are  larger  than  predicted  by  the 
theory  and  also  larger  than  the  values  measured  ca  'he  lover  surface.  As  will  be  discussed  later 
v  chapter  ;  the  differences  observed  have  to  r*»  attributed  to  the  combined  effect  of  airfoil  thickness 
(plus  incidence'  and  the  hcundary  layer.  The  firs*  effect  dominates  on  the  front  part  of  the  airfoil, 
voile  the  boundary  layer  effect  is  .xrn  pronounced  on  the  rear  part. 

\  comparison  between  the  unsteady  pressures  measure-:  cr.  the  upper  surface  cf  the  oscillating 
airfoil  and  the  corresponding  results  : f  'bin  airfcil  theory  frig,  o)  shows  similar  differences  as 
observed  in  quasi -steady  flow.  In  general  the  agreement  r et  veer.  theory  and  experiment  is  reasonable, 
further  it  can  he  noted  that  *  here  is  a  very  satisfactory  agreement  between  the  unsteady  pressures 
measured  directly  with  the  in  situ  transducers  and  tr.e  pressures  obtained  via  the  pressure  tubes. 

F ,  2 . 2  Transonic  flow  with  shock  (  condition  B) 

The  next  example  concerns  oscillations  cf  the  airfoil  about  the  off-design  condition  5  (rig.  u). 
In  this  condition,  being  typical  for  "class! cal”  transonic  flew,  the  upper  surface  carries  a  super¬ 
sonic  region  extending  to  about  50  per  cent  of  the  chord,  which  is  terminated  by  a  relatively  strong 
shock  wave.  As  shown  in  figure  T.I  a  change  ir.  incidence  of  1  degree  results  in  a  sr.ift  of  the  steady 
shock  position  cf  about  1C  per  cent  cf  the  chord.  The  flow  along  the  lower  surface  remains  suberi- 
tical. 

From  the  corresponding  quasi-steady  pressure  distributions  (Fig.  7.II)  it  can  be  deduced  that  along 
the  upper  surface  the  pressure  is  dominated  by  the  effect  cf  the  shock  displacement,  generating  a 
high  pressure  peak,  which  of  course  cannot  be  predicted  by  thin  airfoil  theory.  The  quasi-steady 
pressure  distribution  on  the  subsonic  lower  surface  is  predicted  reasonably  well. 

Unsteady  pressure  distributions  or.  the  upper  surface  are  presented  in  figure  3  for  three  di f fere.t 
frequencies.  These  results  also  show  tr.e  d emir. a r.t  effect  of  the  pressure  peak  due  to  the  moving  shock 
wave.  It  is  noted  that  this  pressure  reax  shifts  from  tr.e  real  part  of  the  pressure  distribution  to 
the  imaginary  part,  with  increasing  frequency.  This  i:;  the  result  cf  the  increased  phase  lag  of  the 
periodical  shock  motion  relative  to  the  metier,  of  the  airfoil,  a  phenomenon  to  be  discussed  in  more 
detail  in  chapter  3.*.i. 

By  representing  the  unsteady  pressure  distributions  in  terms  cf  magnitude  and  phase  angle 
(Fig.  9)  it  further  can  be  shewn  that  the  width  ar.a  the  height  of  the  pressure  peax  associated  with 
the  periodical  motion  of  the  shoe/,  wave  decreases  as  the  frequency  is  increased.  This  is  caused  by 

the  uecrease  cf  the  amplitude  of  the  shock  motion  with  increasing  frequency  (see  also  chapter  3.U.1). 

Concerning  the  phase  curves  in  figure  '*  it  should  be  noted  that  the  measurements  shew  a  jump  of  about 

130  degrees  just  downstream  of  the  mean  position  cf  the  shock  wave.  This  jump  is  present  already  in 

quasi-steady  flow  and  thus  is  not  a  dynamic  effect. 

From  tne  comparison  of  the  measured  pressure  distr i butions  with  the  distributions  calculated  with 
thin  airfcil  theory  (Figs.  3  and  ?)  it  is  evident  tnat ,  as  far  as  the  upper  surface  is  concerned,  this 
theory  is  net  applicable. 


'  i'  special  interest  is  the  unsteady  cenaviour  of  the  airfoil  near  its  "shock-free"  design  con- 
iiricn  (condition  C  in  figure  As  shewn  ir.  figure  I..'. I  a  variation  in  incidence  cf  0.5  degree 
about,  the  design  point  leads  to  a  cons  Herat:  ie  change  of  the  steady  pressure  distribution  along  the 
upr,'*r  surface.  In  particular  ir.  the  supersonic  region,  ranging  from  about  3  per  cent  to  about  t5 
per  e-.r::  of  tne  chord  the  snaps  cf  the  pressure  distribution  changes  considerably.  Further  away  from 
the  design  condition  a  . we ax ;  shoes  wave  shews  up.  At  the  lower  surface  the  steady  pressure  distri- 
:  u*  i on  cr.anges  .ess  drastically.  Notevcrtr.y  is  that  along  the  lower  surface  the  velocity  becomes 
clightiy  our^rcr: ticai ,  tut  still  without  shock  formation. 

The  cr.ang-rs  ir;  steady  distr: cut i or.  result  ir.  a  juasi-steady  distribution  as  given  in  figure  10.11- 
:r.  tne  upper  surface  a  vide  tulgo  -cour  t ,  which  is  caused  by  the  drastic  cr.ange  of  the  pressure  dis- 
*  rir  ution  in  t:.c  cunersor.ic  region.  Most  probably  this  wide  bulge  is  a  feature  ryp.ca-  of  this  *ypv  of 
"chocx-f re^"  airfoil,  characterized  by  a  relative  blunt  nose  and  an  extensive  region  of  supersonic 
fl:v.  A  compar: z-.r.  between  the  measured  quasi-steady  distribution  and  the  curve  determined  with  tnin 
iirfeii  theory  -hows  that  the  : reliction  for  the  upper  surface  is  quite  useless.  For  the  lower  side, 
wr.ere  the  fl  v  remains  almost  succritica. ,  the  differences  between  theory  ar.a  experiment  is  consider- 


A  series  of  fully  unsteady  pressure  distributions  along  the  upper  surface  ir.  terms  of  magnitude 
ar.  i  phase  ar.git?  is  given  in  figure  11.  one  easily  recognises  in  the  magnitude  curves  the  large  r  on  - 
tritutifns  associated  with  the  cr.ar.ges  in  tr.e  snape  cf  the  pressure  distribution  ir.  the  supersonic 
regi  :n  or.  the  front  part  cf  tne  airfcil .  In  addition  a  small  peak  occurs  at  about  c5  per  cent  of  the 
rr.oru,  caused  ry  the  periodical  formation  of  a  weak  shock  in  this  region  see  figure  12).  This  peak 


grtvs  .arger  with  increasing  trequency  as  a  result 


the  increased  strength  of  the  shock  wave.  At 


tne  same  *  im-r  the  bu.*.ge  on  the  front  part  iecreases  with  frequency  and  the  unsteady  pressure  distri¬ 
bution  shown  a  tendency  to  change  ir.  a  direction  towards  the  pressure  distributions  found  for  flow 
condition  3. 

The  phase  curves  shewn  in  figure  .1  behave  very  regular  up  to  about  60  per  cent  of  the  chord.  Then  a 
jump  in  pr.ase  angle  of  about  1:3  degrees  occurs,  which  can  be  attributed  to  the  presence  of  the  shock 


Finally  a  comparison  of  the  measured  unsteady  pressure  distributions  with  thin  airfoil  theory  confirms 
what  could  ie  concluded  already  u;  the  basis  of  the  ]uasi-stcady  data;  for  these  types  of  mixed  flow 
one  :.as  to  rely  . n  other  prediction  methods. 


to  r 
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3.3  Unsteady  aerodynamic  derivatives 

Of  prime  concern  tc  the  aeroelastieia.n  of  course  are  the  ovt-raii  unsteady  aerodynamic  airloads. 
For  this  reason  the  unsteady  aerodynamic  coefficients,  obtained  :y  cn.  r  twine  ir/.v^rat;-.::  f  the 
measured  unsteady  pressure  distributions,  have  been  collected  in  figures  13-1‘.  ,  fv.r  tr.e  characteris¬ 
tic  flow  conditions  A,  B  and  U,  respectively.  For  reference  purposes  l.-.e  results  according  tc  thin 
airfoil  theory  are  given  as  well. 

The  agreement  between  the  theoretical  and  experimental  pressure  distributions  for  the  subsonic 
flow  condition  A  (see  chapter  3.2.1)  is  reflected  also  in  the  curves  cf  figure  13,  representing  the 
unsteady  aerodynamic  derivatives  as  a  function  of  reduced  frequency.  The  largest  deviations,  occuring 
in  the  real  part  of  both  the  normal  force  and  the  moment  derivative,  can  be  attributed  to  the  differ¬ 
ences  in  the  pressure  distributions,  which  do  exist  already  in  quasi-steady  flow  (Fig.  5)- 

For  the  transonic  flow  condition  B  (see  figure  iu)  the  differences  between  tnecry  and  experiment 
are  considerably  larger  than  in  the  preceding  fully  subsonic  example.  This  is  true  also  for  the  un¬ 
steady  derivatives  in  the  "shock-free"  design  condition  C  (Fig.  15).  A  comparison  between  figures 
1U  and  15  learns  that  for  the  design  condition  the  deviations  from  thin  airfoil  theory  are  of  the 
same  order  of  magnitude  as  for  the  "classical”  transonic  flow  condition  B. 

The  behaviour  of  the  aerodynamic  coefficients  in  a  transonic  flow  with  shock  wave  can  be  corre¬ 
lated  qualitatively  with  the  presence  of  the  dominant  pressure  peax  generated  by  the  oscillation  cf 
the  shock.  As  indicated  schematically  in  figure  16  (representing  for  instance  the  results  of  flow 
condition  B),  the  pressure  peax  associated  with  the  shock  wave  is  responsible  for  a  shift  in  unsteady 
lift  and  moment  indicated  by  a  1.  At  small  reduced  frequencies  the  real  part  of  the  normal  force 
derivative,  ka ,  is  larger  than  predicted  cy  theory.  As  the  frequency  increases,  the  real  part  de¬ 
creases  faster  than  the  curve  for  thin  airfoil  Vr.-_-.ry,  while  the  imaginary  part  becomes  suer,  more 
negative  than  predicted.  This  behaviour  10  c  .-re.atei  with  tr.e  shift  of  the  pressure  peax  due  to  the 
shock  from  the  real  part  to  the  imaginary  port  cf  the  -unsteady  pressure  distribution,  as  r.as  been 
shown  in  figure  0  (see  also  section  3.2.-'. 

The  same  phenomenon  is  responsible  for  the  change  indicated  by  1,  in  the  moment  derivative,  .  The 
remaining  part  of  the  deviation  in  the  moment  derivative  is  caused  by  the  circunstance  that  the 
mentioned  pressure  peak  is  located  downstream  of  the  quarter  chord  point,  thus  giving  rise  to  a  rear¬ 
ward  shift  of  the  aerodynamic  centre.  For  the  present  example  this  shift,  expressed  as 


c  2  k 

a 


can  be  estimated  roughly  at  5  per  cent  of  the  chord. 

In  the  figures  13-15  results  are  given  for  the  airfoil  with  and  without  transition  strip.  For 
flow  conditions  A  and  B  no  significant  difference  is  observed.  However,  in  the  delicate  "shock-free" 
design  condition  C  the  flow  is  more  susceptible  to  disturbances  over  the  front  part  of  the  airfoil 
and  thus  more  sensitive  to  the  presence  of  the  strip.  For  a  more  detailed  account  on  this  sensivity 
reference  is  made  to  reference  T. 

Further  it  should  be  remarked  here  that  the  measured  data  are  given  without  tunnel  wall  correct¬ 
ion,  since  reliable  methods  to  determine  this  effect  in  unsteady  wind  tunnel  tests  are  not  yet 
available.  An  estimate  of  the  amount  of  wall  interference  involved  in  the  present  tests  will  be  given 
in  chapter  it  on  the  basis  of  some  quasi-steady  flow  calculations. 

3.1*  Remarks  on  the  unsteady  shock  wave  motion 
3.l*.l  Effect  of  frequency 

With  the  help  of  optical  flow  studies  additional  information  is  obtained  about  the  periodical 
motion  of  the  shock  waves  in  flow  condition  B.  From  figure  17,  giving  the  time  histories  of  the  shock 
displacement  for  different  frequencies,  it  follows  that  the  shock  wave  performs  nearly  sinusoidal 
motions  (similar  to  the  type  A  motion  described  in  reference  8).  Further  the  phase  lag  of  the  shock 
motion  relative  to  the  airfoil  motion  increases  with  frequency,  while  the  amplitude  of  the  shock 
motion  decreases.  The  latter  corresponds  very  well  with  the  observations  mentioned  earlier  concerning 
the  contribution  of  the  moving  shock  wave  to  the  unsteady  pressure  distributions  (see  Figs.  8  and  9). 

A  closer  examination  of  the  phase  lag  of  the  shock  motion  with  respect  to  the  airfoil  motion 
(Fig.  18)  learns  that  an  almost  linear  relationship  exists  between  frequency  and  phase  lag.  This  im¬ 
plies  that  there  is  a  constant  time  lag  between  the  motion  of  the  airfoil  and  the  shock  wave  motion. 
In  relation  to  this  it  is  of  interest  to  recall  the  investigation  of  Erickson  and  Stephenson  (Ref.  9) 
who  have  found  that  a  fixed  relation  seems  to  exist  between  the  phase  lag  of  the  shock  motion  and  the 
time  required  for  a  pressure  impulse  to  travel  from  the  trailing  edge  to  the  shock  wave.  Indeed  this 
travelling  time  seems  to  be  a  logical  parameter  for  an  airfoil  with  a  large  supersonic  region  ter¬ 
minated  by  a  shock  wave,  because  this  is  the  time  period  after  which  major  changes  in  flow  condition, 
namely  changes  in  flow  direction  at  the  trailing  edge  (Kutta  condition)  can  be  felt  by  the  shock  wave 
(Fig.  19). 

The  time  required  to  forward  information  from  the  trailing  edge  to  the  shock  wave  amounts 
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with  Mqoc  being  the  local  Mach  number  and  aloc  the  local  velocity  of  sound.  Tue  to  the  gradient  ir. 
Mach  number  normal  to  the  airfoil  surface  the  acoustic  waves  propagate  along  paths  away  from  the 
airfoil.  Therefore  the  propagation  speed  in  upstream  direction  will  be  some  average  between  the  value 
of  (l-Mloc)  aloc  near  the  airfoil  surface  and  the  free  stream  value.  To  account  for  this  effect  the 
following  value  of  the  local  Mach  number  has  been  introduced: 


Mloc*  P  |Mloc(at  the  surface)-  Mx|+ 


v i t H  Veing  a  relaxation  factor ,  which  \.'is  a  value  :  ■nvr...n  0  a: .a  1.0. 

V  > r  an  excellent  agreement  :u  :  etveon  *  ho  travel  time  of  *.*:<■»  ”Y. uf-.a  vav**:;"  and  tne 

e  .:rre?p,:ndir;g  phase  lag  *•:*  •  r. 'tur. ,  as  i -.iemor.r  t  rati- i  : :'.gur“  i.-;.  "ml;  valae  f  r  agree.: 

very  veil  with  the  relaxation  factor  ap:l:*-.i  i  r.  •/>■'  LoufciM  ;  aft  ice  .retr.oj  of  L’LR  -**e 

Hoes,  Kef.  1C ) ,  in  which  ir.  a  remi-emnir i Mil  vav  a  local  Much  number  :errectitr.  Is  ; r.trt -.iuc«-d .  I:. 
tr.is  method  also  the  factor  R  is  used  to  account  for  tr.e  gradients  in  Maci:  number  normal  to  the  air¬ 
foil  surface. 
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wave  it  has  leer,  cc served  that  the  shock  me 
of  the  shock  metier,  is  almost  prepert ior.al 
\  see  for  instance  figure  LI ) .  This  makes  it 


ion  taxes  place  almost  sinusoidal  and  that  the  amplitude 
c  the  amplitude  of  the  sinusoidal  motion  of  the  airfoil 
possible  to  introduce  the  schematized  moael  cf  figure  21 
in  which  the  change  in  pressure  in  a  fixed  point  A  is  considered,  while  the  shock  wave  performs  a 
sinusoidal  motion  of  amplitude  m. 

As  derived  in  reference  6  the  local  shock  strength  in  point  A,  located  within  +  he  sheen  trajectory, 
can  be  written  as: 
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denotes  the  unit  step  function  and  dp  the  variation  in  shock  strengnt  during  the 


shock  wave  notion.  For  strong  shock  waves  and  small  amplitude  motions  the  last  term  ir.  the  above 


expression  car.  re  discarded  relative  to  (p-.-p.)  .  When  ( po— Pc  )  x  '*s  described  as  a  function  of  time 


a  block  type  signal  occurs  (see  figure  22),  cf  whicr.  the  Fourier  decomposition  yields: 
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The  corresponding  distribution  of  the  first  four  Fourier  components  along  the  trajectory  of  the  sr.ee 
action  are  shown  ir.  figure  23.  From  the  distribution  of  the  mean  steady  value  it  follows  that,  due  t 
‘he  oscillatory  metier,  of  the  shock  wave,  the  jump  i..  the  steady  pressure  distribution  is  spread  cut 
ever  the  shock  trajectory.  The  distribution  of  the  component  with  the  same  frequency  as  the  airfoil 
metier,  chews  a  maximum  of  2/t  times  the  steady  pressure  jump  (pj-pq) 

x? 

Integration  cf  the  various  components  over  the  shock  trajectory  tc  obtain  the  contribution  to  the 
overall  unsteady  lift  and  moment  learns  that  the  lift  contains  only  a  contribution  of  the  fundamental 

frequency.  The  resulting  unsteady  moment  also  contains  a  term  (p,-p.  )  .  .  cos  2<ut .  So  it  can  Ve 

*-  o 

expected  that  the-  second  harmonic  shews  up  first  ir.  the  unsteady  moment. 

From  the  considerations  given  above  it  fellows  also  that  measuring  the  first  Fourier  component 
of  the  pressure  signals,  as  is  acne  in  the  present  tests  via  the  tubing  system,  gives  by  chordw-.se 
integration  a  correct  value  of  the  unsteady  lift.  As  far  as  the  moment  is  concerned  t.ue  second  har¬ 
monic  of  order  x~  can  r.ot  Ve  distinguished. 

EXFECTZD  CAPABILITY  CF  THE  HEW  CALCULATI oil  METHODS 

Ir.  the  preceding  analysis  cf  the  experimental  data  the  results  of  thin  airfoil  theory  have  teen 
added  as  a  reference  for  two  reasons.  Firstly  these  results  serve  as  a  simple  basis  for  the  distinct 
of  tr.e  typical  transonic  phenomena  and  secondly,  linear  lifting  surface  theory  is  widely  used  ir.  aer 
elastic  applications.  As  long  as  the  flow  is  moderately  subsonic  thin  airfoil  theory  has  proven  to 
be  a  rather  adequate  -col  indeed.  However,  from  tr.e  preceding  discussions  it  is  apparent,  that  cal¬ 
culation  methods  for  transonic  flow  should  include  the  effect  of  airfoil  thickness,  incidence  and  - 
if  shock  waves  are  present  -  also  the  effects  of  'he  periodical  shock  wave  motion. 

In  recent  years  considerable  progrecs  has  beer,  achieved  in  solving  the  non  linear  equa.'ons  for 
unsteady  transonic  flow  (reviews  on  the  eurrer.‘  status  are  giver,  in  references  12-20). 

With  one  exception  (Ref.  21)  all  new  calculation  methods  are  dealing  with  inviscid  flow.  In  order  to 
get  an  impression  about  the  improvements  one  might  expect  from  these  methods  some  comparisons  will  b 
presented  between  theory  and  experiment  for  the  NLR  7301  airfoil.  At  this  moment  the  comparison  is 
limited  to  quasi-steadv  flow,  but  in  the  near  future  comparative  studies  will  be  performed  also  for 
fully  unsteady  flow. 

Considering  first  the  quasi-steady  cane  has  the  advantage  that  a  reasonable  estimate  can  be 
given  of  the  effect  of  the  boundary  layer,  by  using  an  existing  method  for  steady  transonic  flow, 
which  includes  the  displacement  effect  of  the  boundary  layer  (Bauer,  Korn,  Garabedian  and  Jameson 
(Ref.  22)).  Further  for  quasi-steady  flew  a  rather  accurate  estimate  can  be  given  of  the  severity 
of  interference  from  the  slotted  tunnel  walls,  a  situation  which  is  not  yet  reached  for  unsteady 
measurements . 


Mach  number  as  in  the  experiments  ,  rut  for  un  i  no  wen':*-  u.,  ..  -  j  j«'  :•  *  ■  rr'-:"  ;  ,r  wuj  . 

interference.  For  the  ULR  Pilot  tunnel  thio  correct i nas  t  e**n  •  •j*.  :i;.  .  ' 

Ja=  ;.  ( 

where  C_  is  steady  lift  coefficient  and  C  a  coefficient,  depending  the  frv»?  stream  M-aor.  r.umtor 
(Fig.  2$. 

It.l  Fully  subsonic  Flow  (condition  A) 

Calculated  and  measured  results  for  the  airfoil  it  the  subset:;  flow  condition  A  are  sr.owr.  ir. 
figure  25.  Figure  25.1  reveals  a  significant  effect  of  the  counuary  layer  in  the  steady  tear,  posi¬ 
tion  of  the  airfoil.  The  corresponding  quasi-steady  results  ; Fig.  15.11)  demonstrate  • nat  the  devia¬ 
tions  of  the  test  results  from  thin  airfoil  theory,  as  discussed  ir.  chapter  b.2.1,  are  uue  to  the 
combined  effects  of  thickness,  incidence  and  viscosity.  The  effect  f  thickness  and  incidence 
dominates  on  the  front  part  of  the  airfoil  and  the  effect  cf  viscosity  towards  the  rear. 

For  the  quasi -steady  results  the  wall  correction  has  seen  applied  on  the  measured  data,  since 

this  effect  can  be  translated  simply  in  an  additional  change  of  effective  incidence  due  to  t  r.e 
change  in  lift.  This  additional  change  has  to  be  substracted  from  the  geometrical  change  in  incidence. 
From  the  results  in  figure  -.11  ir  ear  be  noted  that  the  tunnel  walls  r.ave  a  considerable  effect. 

*4. 2  Transonic  flow  with  shock  wave  (condition  3) 

The  second  example  deals  with  the  transonic  flew  condition  5  (Fig.  Ir.  steady  flow  (Fig.  26.1) 
viscosity  again  has  a  large  effect,  in  particular  or.  the  location  of  tne  shock  wave.  The  importance 
of  inserting  boundary-  layer  effects  is  reflected  also  ir.  the  quasi-steady  results  trig.  26.11).  In 
the  upper  surface  a  considerably  improved  prediction  is  obtained,  when  thickness  and  boundary  layer 
effects  are  considered  simultaneously.  Especially  the  location  cf  the  high  pressure  peak  resulting 
from  the  shift  in  shock  position  is  predicted  much  better. 

The  improvements  achieved  cam  be  observed  also  in  the  quasi-steady  aerodynamic  coefficients, 
collected  in  table  1.  For  instance  these  data  show  that  a  ”x=  0.7  thickness  ar.d  incidence  are  res¬ 
ponsible  for  an  increase  of  the  thin  airfoil  value  of  the  normal  force  coefficient,  ka ,  of  mere  than 

50  per  cent.  The  inclusion  of  the  boundary-layer  leads  to  a  decrease  of  the  order  of  35  per  cent, 
as  can  be  observed  by  comparing  the  results  with  and  without  boundary  layer,  both  obtained  with  the 
non-conservative  calculation  scheme  (the  conservative  scheme,  which  guarantees  the  best  numerical 
solution  of  the  transonic  flow  equations  did  not  converge  for  inviscid  flow  so  this  value  could  r.ct 
be  added).  From  the  last  two  columns  it  follows  that  the  tunnel  wall  effect  ir.  tne  present  tests  is 
considerable  and  accounts  to  about  25  per  cent.  At  Moc=  0.5  the  effects  mentioned  are  less  than  at 
transonic  speed,  but  still  significant. 


TABLE  1 

Quasi-steady  aerodynamic  derivatives  (MLR  7301  airfoil) 


Thin 

airfoil 

theory 

Inviscid  theory 
+  thickness 

Inviscid  theory 
+  thickness 
+  boundary  layer 

Inviscid  theory 
+  thickness 
+  boundary  layer 
+  wall  interference 

Experiment 

Non¬ 

conservative 

F-D  scheme 

Conservative 

F-D  scheme 

a 

k 

m 

K. 

m 

m 

m 

k 

m 

k 

m 

0 

a 

a 

a 

a 

a 

a 

a 

a 

a 

0.5 

0.85° 

2.31 

3 

2.73 

1 

2.53 

-0.036 

2.22 

-0.032 

2.18 

-0.090 

0.7 

3.00° 

2. 8o 

0 

9. 29 

0.11 

_ 

3.21 

0.00 

3.92 

_ 

-0.22 

_ 

3.23 

-0.18 

3.20 

-0.39 

F-D=  Finite  difference 


From  the  examples  discussed  so  far  a  good  impression  is  obtained  about  the  improvements  which 
can  be  expected  at  most  from  the  inclusion  of  thickness  and  incidence  theories.  Clearly  the  inclu¬ 
sion  of  these  effects  is  an  important  step  forward,  which  on  itself,  however,  does  not  lead  to  im¬ 
proved  predictions.  A  genuine  improvement  in  this  respect  can  be  achieved  only  if  the  second  step  is 
made  also,  i.e.  the  inclusion  of  boundary  layer  effects. 

A  weak  point  in  the  considerations  given  above  is  seemingly  that  the  examples  deal  with  a  rela¬ 
tively  low  Reynolds  number  (  -2.10^).  However,  similar  calculations  for  higher  values  of  this  param¬ 
eter  (up  to  30.10^,  with  fixed  transition  point)  do  not  exhibit  a  significant  sensivity  to  Reynolds 
number  changes.  This  seems  to  indicate  that  under  full  scale  conditions  the  effect  cf  viscosity  re¬ 
mains  of  the  same  order  of  magnitude  as  shown  here. 

U . 3  The  "shock- free"  design  condition  (condition  C) 

To  conclude  the  evaluation  of  the  capability  of  advanced  theories  cn  the  basis  of  quasi-steady 
flow  the  "shock-free  flow  condition  C  will  be  considered.  For  this  purpose  a  comparison  is  made  be¬ 
tween  results  calculated  for  the  theoretical  "shock-free”  design  condition  ana  results  measured  for 
condition  at  which  "shock-free"  flow  is  obtained  ir.  the  wind  tunnel. 

In  this  way  the  circumstance  that  the  experimental  design  condition  (i.e.  Mach  number  and  incidence) 
differs  from  the  inviscid  theoretical  design  condition  can  be  discarded,  assuring  that  both  theory 
and  experiment  deal  with  the  carefully  balanced  condition  of  "shock-free"  flow. 

The  steady  pressure  distributions  computed  for  incidences  at  and  around  the  design  condition 
(rig.  27.1)  exhibit  in  the  supersonic  region  at  the  upper  surface  the  same  ir.ar.ied  changes  in  the 
share  of  the  pressure  distribution  as  observed  in  the  measurements  'Fig.  10).  The  lower  surface  be- 


thin  airfoil  theory.  In*-1  r 
dieted  reasonably  v*-*  L 1  uno 
e xpe o  t  at i on  that  met roa c.  : 
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From  the  preceding  vvaluat  i  ,n  it  in  apparent  that  ?  r>*  inclusion  f  airh'ii  *  hi cknvss ,  ir.ciuer.ee 
ana  transonic  shoe*  motions  in  inviscid  flow  calculations  leads  to  an  improvement  :f  tre  theoretical 
predictions  in  an  at  i_is’  juacitative  cer.oe.  in  quantitative  sense  i  .urge  discrepancy  vitn  The  real 
flow  will  remain  as  a  result  of  the  boundary  layer,  vhicn  tc  a  ..urge  extent  determines  tne  final  lo¬ 
cation  cf  the  shock  and  by  that  the  overall  unsteady  airliads.  however,  as  the  modelling  cf  unsteady 
boundary  layers  is  only  in  its  first  phase  for  a  review  cf  *  he  present  status  reference  is  made  to 
Ref.  1;  it  is  ur.lix- ly  that  i:.  the  near  future  s. ;  crust  i  sat  ou  cal  ruiat  i  on  methods  will  receme 
available  for  this  purpose.  Iherefcre  in  tne  coming  penou  an  engineering  type  of  v.  preach  has  to  he 
f:il_v-u.  In  this  respect  t:.e  ideas  developed  by  Magnus  ar.G  7 ;  sn  mara  Re:'.  1 ;  or- serve  attention, 
since  ri.  ir  relatively  simple  "viscous  ramp"  model  lends  itself  for  easy  : mr  1  oment at : on  in  ir.viscii 
cal cul at i or.  met h.oas . 

From  a  eonputat icnal  point  if  view  small  perturbation  methods  are  very  attractive.  It  should 
be  investigated,  therefore,  what,  the  limits  if  such  methods  are,  in  parti  .-u^ur  wnern  app.ie-i  to  t  h  i;  s 
supercritical  airfoils  of  the  type  as  the  considered  in  this  paper.  In  *nis  respect  it  should  be 
noted  that  the  impressions  given  in  the  preceding  eaapttrr  about  the  improvements  attainable  with  the 
new  calculation  methods  are  based  on  solutions  cf  the  full  potential  equation,  without  assuming  small 
perturbations . 

Further  from  the  considerations  cf  the  quasi- steady  results  of  the  N  LR  Toll  airfoil  it  has  be¬ 
come  clear  tnat  in  order  to  improve  the  reliability  of  comparisons  between  theory  and  wind  tunnel 
data  there  is  an  urgent  need  for  methods  to  assess  the  amount  of  wail  interference  in  unsteady  ex¬ 
periments  in  transonic  test  sections  with  slotted  or  pereus  walls.  Finally  the  insight  with  respect 
to  the  eff--o*  of  Reynolds  number,  being  already  a  crucial  parameter  in  steady  transonic  flew,  should 
be  increased  by  performing  tests  in  a  high  Reynolds  number  test  facility. 
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I  O  UPPER  SURFACE  - 


Fig.  7  Steady  and  quasi  -  steady  pressure  di stributions  on 
NLR  7301  airfoii  in  transonic  flow  with  shock  wave 
(Condition  B ) 


Fiq.  8  Deve/oo'nenf  of  unsteady  orcssure  distribution*  wth 
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Fig.  9  Effect  of  a  shock  wove  on  the  unsteady  pressure 
di  stnbut  ion  s  (Condifion  B) 
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Fig.  13  Unsteady  normal  force  ond  moment  coefficients  os  o 
function  >f  frequency  m  subsonic  flow  l  Condition 


Fig,  14  Unsteady  normal  force  and  moment  coefficients  os  o 
function  of  frequency  (Condition  8) 
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F<q.  17  Periodical  shock  wen  nit  o  >  <.-)  dependence  of 
f'emency  <Cond  tion 


Fig.  26  -  I  Steady  pressure  d>  stributions  <■'  mean  position 
‘Condition  8/ 
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Fig.  27  -II  Quasi  -  steady  pressure  distributions  for  the  theoretical 
"  shock  -  free"  design  condition  fCondif«on  C) 


